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Integral  transform  formulae  are  given  for  the  time-harmonic  acoustic' 
pressure  and  particle  velocities  from  which  the  acoustic  intensity  vectors 
in  the  interior  and  exterior  fluid  may  be  calculated  numerically.  Some 
plots  of  interior  intensity  vectors  due  to  point-force  or  point-source 
excitation  demonstrate  their  potential  usefulness  as  a  visual  aid  to  the 
understanding  of  wave  propagation  in  fluid-filled  pipes. 
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LIST  OF  SYMBOLS 

Cartesian  and  cylindrical  coordinates 

axial,  tangential  and  radial  displacements 
at  shell  mid-surface 

axial  and  circumferential  harmonic  wavenumbers 

Fourier  transforms  of  shell  displacements 

interior  and  exterior  pressures 

interior  acoustic  particle  displacements 

exterior  acoustic  particle  displacements 

acoustic  intensity 

radian  frequency  of  vibration 

density  and  sound  velocity  of  interior  fluid 

density  and  sound  velocity  of  exterior  fluid 

wavenumber  in  fluid,  w/c^ 

radius,  thickness  and  density  of  shell 

=Eh/(l-a  ).  E  is  Young's  modulus 
c7  is  Poisson's  ratio 

=h2/l2a2 

=1  when  n=0,  and  2  otherwise 

amplitude  of  point  force  located  at 
cylindrical  coordinates  (a, 0,0) 

amplitude  of  point  source  located  at 
cylindrical  coordinates  (x^,0, 0) 

2  2  p  p  p  2 

=  (x-xQ)  +y~+z  =  r  +x^-2rx3cos(^)+z 

3essel  functions  and  their  derivatives, 

H=J+iY 
n  n  n 

2.  2 

(k^-at  )2  with  lm.(b^)>0  to  satisfy  the 
radiation  conditions 

hyseretic  loss  factors  of  shell  and  fluid 


velocity  and  complex  conjugate  of  pressure 


INTRODUCTION 


Some  recent  work  by  Fuller  and  Fahy  Llj  has  contributed  to  the 
understanding  of  the  physics  of  wave  propagation  in  piping  systems 
through  their  discussions  on  the  physical  interpretation  of  the  disp¬ 
ersion  relation  and  the  energy  distribution  of  free-waves  between 
fluid  and  shell.  James  [2, 3 J  has  shown  that  the  dispersion  relation, 
presented  as  wavenumber  versus  frequency  plots,  helps  to  promote  a 
qualitative  understanding  of  sound  radiation,  power  radiation  and 
near-field  responses. 

Publications  by  Fahy  [4,5]  have  stimulated  interest  in  the  meas¬ 
urement  of  acoustic  intensity  using  the  two-microphone  technique.  The 
paper  by  Reinhart  and  Crocker  [6]  provides  an  instructive  account  of  a 
practical  application.  The  usefulness  of  intensity  vector  plots 
obtained  from  theoretical  work  has  been  demonstrated  by  Spicer  L7-  for 
the  particular  case  of  an  elastic  plate  with  mass-spring  attachments 
undergoing  forced  vibrations.  His  plots  illustrate  neatly  the  physics 
of  the  fluid-plate  interaction  in  terms  of  energy  flow  in  the  fluid;  in 
particular,  the  fluid  subsonic  surface  wave  and  angled  leaky -wave  beam 
are  clearly  evident. 

The  main  purpose  of  the  work  contained  herein  is  to  discuss  a 
limited  number  of  wave-frequency  and  interior  acoustic  intensity  vector 
plots  for  the  case  of  an  infinite  cylindrical  elastic  shell  under  the 
action  of  point-force  or  interior  point-source  excitation.  It  is  hoped 
that  intensity  vector  plots  will  display  features  that  enhance  the 
physical  interpretation  of  the  wavenumber-frequency  plots.  The  plots 
may  also  be  of  interest  to  experimenters  who  are  using  the  technique  of 
acoustic  intensity  measurement  to  investigate  sound  transmission  in  pipes 
because  a  combination  of  interacting  experimental  and  theoretical  work 
is  needed  to  achieve  noise  control  in  complex  pipe  networks. 

2.  PROBLEM  FORMULATION 

(a)  General 

An  infinite  cylindrical  thin-walled  shell  contains  and  is  immersed 
in  an  acoustic  fluid.  The_displacements  of  the  shell  wall  satisfy  an 
eighth-order  shell  theory  LS],  and  the  interior  and  exterior  sound,  press¬ 
ures  satisfy  acoustic  wave  equations.  The  time-harmonic  excitation  is 
provided  by  mechanical  point-forces  at  the  shell  surface  or  acoustic 
point-sources  in  the  interior  fluid.  The  time  factor  exp(-iwt)  is 
omitted  from  all  equations.  The  geometry  is  shown  in  Figure  1.  The 
mathematics  is  presented  with  a  minimum  of  elaboration  and  explanation 
because  much  of  it  has  been  given  in  greater  detail  elsewhere  [2, 33. 


(b)  Pressure  and  Particle  Displacements  Inside  Pipe 


The  excitation  is  assumed  to  be  of  even  ^-dependence  so  the 
pressure,  radial  displacement  and  axial  displacement  can  be  expressed 
as  the  Fourier  transforms 
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The  spectral  pressure  is  of  the  form  [2, 3D 


p^(r,n,0t)  =  /,1w2[jn(b1r)/b1Jl(b1a)]w(n,a)  +  A(r,n,a)  (2.2) 


where  W(n,ot)  is  the  radial  spectral  displacement  at  the  shell  surface  and 
A(r,n,ot)  is  identically  zero  for  point-force  excitation.  A  matrix  formula 
for  ft(n,0t)  for  the  special  case  of  Goldenveizer-Novozhilov  shell  theory 
[3D  is  given  in  the  Appendix.  For  point  source  excitation  [2, 3D 


S(r,n,a)  =  ^Poen[jri(b1xo)/j;(b1a)Ijn(b1r)Y;(b1a)-Yn(b1r)j;(b1a)Df  r>xQ 

(2.3) 

"  ^oenCjilblr)/jn(bla)ljn(Vo^n(bla)‘Yn(blxO)Jn(bla)]'  r<xo 

The  radial  and  axial  acoustic  particle  displacements  are  related  to 
the  pressure  by  the  formula 


[dPj/dr,  dPx/dzD  =  P1w2[ur]_,  uzlD  (2 .'+) 

from  which  the  spectral  displacements 

uri(r,n,a)  =  ^Jn(bir)/Jn(b1a)]w(n,a)  +  ( 0Lw2)"1dS(r,n,a)/dr  (2 .5) 

uzl(r,n,a)  »  iaCjii(b1r)/b1J^(b1a)Dw(n,a)  +  ia(  ^w2)"1!^,  n,a)  (2.6) 


may  be  obtained  by  the  use  of  equations  (2.I-2.3).  In  equation  (2.5) 

*VdT  *  T^oe/Jn(blxo)/JA(bla)^JA(blr)Yn(bla)*YA(blr)jA(bla):i»  r>xo 
-  ^0«n[JA(b1r)/j;(b1a)Db1[jn(b1xo)Y;(b1a)-Yn(b1xo)j;(b1a)D,  r<x 


(2.7) 


3.  NUMERICAL  EVALUATION  OF  INTEGRALS 


The  integrals  in  equations  (2.12)  must  be  evaluated  numerically 
because  closed-form  expressions  are  not  available.  Their  numerical 
approximation,  on  truncating  the  infinite  limits  to  finite  values,  is 
based  on  a  simple  adaptive  Gaussian  quadrature  scheme  of  order  2  that 
has  been  used  L3j  to  evaluate  similar  integrals.  Because  the  Fourier 
integral  transform  representation  is  arbitrary  with  respect  to  a 
solution  of  the  homogeneous  equations,  it  is  necessary  to  introduce 
damping  into  the  system  before  numerical  evaluation  is  possible.  The 
usual  procedure  is  adopted  here  to  include  dissipation,  viz  set  2=E(l-iPs) 
and  c-c(l-i'L)  where  n  and  r<  are  the  hysteretic  loss  factors  in  the 
shell  wall  and  fluid  respectively, 

When  the  excitation  is  a  radial  force  or  source  located  at  z=0, 
then  the  computational  times  may  be  halved  by  writing  equations  (2.12) 
in  the  form 
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It  is  also  evident  that  further  savings  in  computer  time  are  possible  by 
parallel  computation  at  an  array  of  2-values,  due  to  the  simplicity  of 
the  z-dependent  term  in  relation  to  the  r-dependent  term.  Methods  based 
on  the  fast  Fourier  transform  algorithm  are  unlikely  to  be  of  much 
greater  efficiency. 

4.  NUMERICAL  EXAMPLES 

(a)  General 

For  the  particular  case  of  a  water-filled  shell  surrounded  by  a 
vacuum.  Figures  2-3  show  wave-number-frequency  plots  and  Figures  4-10 
show  intensity  vector  plots  for  the  n-0  and  2  harmonics  separately.  The 
material  and  geometric  constants  in  31  units  that  have  been  used  in  the 
computations  are  as  follows: 

steel:  3-19. 5xl010  v=0.29  P-7700.0  h-0.01  a-0.10  n  =0.02 

s 

water:  P =1000.0  c=1500.0  ^=0.001 

The  radial  point-force  was  located  at  the  cylindrical  coordinates  (a, 0,0) 
and  the  interior  point-source  was  located  at  the  point  (2a/3,0,0).  The 
fluid  loss-factor,  P^-0.001,  is  sufficiently  large  enough  to  allow 
numerical  integration  of  equations  (3  .1)  without  ill-conditioning,  yet  it 
is  sufficiently  small  to  avoid  significant  attenuation  in  the  fluid  over 
the  maximum  axial  distance  considered,  viz  0.4m. 

The  intensity  vector  plots  are  3hown  for  the  separate  values  of  the 
circumferential  harmonic,  n,  in  order  to  aid  interpretation  via  the  wave- 
number  versus  frequency  plots;  also,  the  individual  intensities  have 


physical  meaning  because  they  are  measureable.  Due  to  the  symmetry  of  the 
excitations,  the  intensity  vectors  are  plotted  for  z>0  only;  plots  for 
2<0  axe  simply  mirror  images.  In  all  Figures  the  plotted  lengths  of  the 
normalised  vectors  are  proportional  to  ). 

(b)  Wavenumber-Frequency  Plots 

The  wavenumbers  for  free-wave  propagation  in  the  axial  direction  are 
those  real  values  of  Cl  for  which  the  determinant  of  the  matrix  equation  (Al'! 
vanishes.  Plots  of  these  values  versus  frequency,  for  a  selected  value  of 
n,  are  called  axial  wavenumber  versus  frequency  plots.  The  physical  inter¬ 
pretation  of  these  plots  is  described  elsewhere  L1J;  in  particular,  the 
charging  nature  of  each  branch  with  frequency  is  of  central  importance  as 
its  energy  ratios  between  shell  and  fluid  may  alter  considerably. 

Figure  2  contains  the  wavenumber-frequency  plots  of  a  steel  pipe 
whose  interior  and  exterior  are  vacuums.  Superimposed  on  these  plots  are 
the  wavenumbers  of  a  rigid-walled  circular  duct  of  water.  Figure  3 
contains  the  plots  for  a  steel  pipe  containing  waters  briefly,  for  n=0 
the  branches  labelled  1-4  cut-on  close  to  a  plane  fluid  wave,  an  axial 
shell  wave,  a  radial  shell  wave  and  a  fluid  wave,  respectively  -  the 
branch  'O’  is  a  wave  of  pure  torsion;  for  n=2  the  branches  labelled  1-3 
cut-on  close  to  a  radial  shell  wave,  a  fluid  wave  and  a  torsional  wave, 
respectively . 

(c)  Rigid  Shells  Source  Excitation 

Figures  4a  and  4b  show  the  intensity  vectors  of  the  n=0  harmonic 
at  frequencies  of  4  and  12kHz,  respectively.  The  plane  wave  mode  alone 
propagates  at  4kHz,  while  at  12kHz  it  is  just  evident  that  the  two  prop¬ 
agating  modes  are  combining  to  form  a  'weak'  interference  pattern.  The 
intensity  vectors  on  the  shell  axis  must  be  parallel  to  that  axis  because 
the  fluid  particle  velocity  in  the  radial  direction  vanishes  there. 

% 

Figures  5a  and  5b  show  the  intensity  vectors  of  the  n=2  harmonic 
at  frequencies  of  4  and  12kHz,  respectively.  At  4kHz  the  intensity 
vectors  decay  rapidly  with  distance  because  of  the  absence  of  a  prop¬ 
agating  mode,  while  at  12kHz  there  is  a  single  propagating  mode.  The 
intensity  vectors  on  the  shell  axis  must  vanish  because  both  pressure 
and  particle  velocities  are  identically  zero  there. 

(d)  Steel  Shell:  Source  Excitation 


Figures  6  and  7  show  the  intensity  vectors  of  the  n=0  harmonic  at 
the  frequencies  of  7.2 5  and  12kHz,  respectively.  A  plot  at  4kHz  is  not 
shown  because  it  has  much  the  same  appearance  as  the  plot  of  Figure  4a, 
except  that  the  vectors  at  the  shell  centre  have  marginally  smaller 
magnitudes.  At  7.25kHz,  which  is  just  above  the  cut-on  frequency  of  the 
radial  shell  mode,  modal  interference  effects  and  areas  of  energy  circ¬ 
ulation  are  clearly  evident;  also,  there  is  an  interchange  of  energy 


between  the  fluid  and  the  shell.  Similar  remarks  apply  to  the  plot  at 
12kHz.  Schultz C9]  notes  that  areas  of  acoustic  energy  circulation  are 
the  rule  rather  than  the  exception  in  acoustic  fields;  they  are  even 
present  in  standing-wave  fields  due  to  the  superposition  of  modes 
whose  individual  intensities  vanish! 

Figure  3  shows  the  intensity  vectors  for  the  n=2  harmonic  at  -kHz 
where  there  is  a  single  propagating  mode.  Most  of  the  energy  from  the 
source  immediately  enters  the  shell,  and  the  propagating  fluid  surface 
wave  associated  with  a  radial  shell  wave  is  evident.  At  12kHz,  not 
shown  here,  the  n=2  plot  is  little  different  from  the  rigid  shell  plot 
of  Figure  5b. 

(e)  Steel  Shell:  Force  Excitation 

Figure  9  is  the  plot  of  the  n=0  harmonic  at  7.25kHz .  It  has  many 
features  in  common  with  Figure  6  (the  plot  with  source  excitation)  with 
the  interchange  of  energy  between  fluid  and  shell  being  particularly 
marked.  Figure  10  is  the  vector  plot  of  the  n=2  harmonic  at  12kHz.  In 
the  range  r<a/2,  a  comparison  with  Figure  5b  shows  that  the  energy  is 
due  to  an  essentially  fluid-type  mode;  however,  in  the  range  r>a/2, 
interference  between  shell  and  fluid  modes  results  in  a  considerable 
interchange  of  energy  between  fluid  and  shell. 

5-  CONCLUDING  REMARKS 

Formulae  have  been  given  from  which  the  acoustic  intensity  vectors 
interior  and  exterior  to  an  elastic  shell  may  be  calculated  numerically. 
Enthusiasts  of  wave  propagation  problems  will  find  the  work  of  Fuller  and 
Fahy  ClJ  a  valuable  aid  to  their  interpretation.  The  limited  number  of 
plots  of  interior  vectors,  due  to  point- force  or  point-source  excitation, 
that  are  presented  are  of  sufficient  interest  to  justify  further  numer¬ 
ical  work  to  enhance  their  physical  interpretation  and  their  practical 
usefulness  in  noise  control. 

First,  it  is  necessary  to  obtain  intensity  vector  plots  for  an 
extended  range  of  parameters  that  includes 

(i)  sufficiently  small  frequency  steps; 

(ii)  the  individual  (n=0,l,2)  harmonics  together  with  their  sum; 

(iii) the  presence  of  an  exterior  fluid; 

(iv)  a  range  of  shell  materials. 

This  extended  range  of  plots  together  with  wavenumber  versus  frequency 
plots  should  assist  the  understanding  of  wave  propagation  in  and  sound 
radiation  from  fluid-filled  pipes. 

Secondly,  intensity  vectors  in  the  presence  of  point  or  axisymmetric 
constraints  on  the  shell  surface  may  be  of  practical  value  to  investig¬ 
ations  into  wave  transmission  phenomena  in  complex  piping  networks  in 


which  non-uniformities  are  common.  No  new  procedure  is  involved  here; 
the  method  of  dynamic  stiffness  coupling  can  be  used  to  determine  the 
reaction  forces  caused  by  the  constraints,  whence  the  determination 
of  the  intensity  vectors  proceeds  in  an  uncomplicated  way.  The  procedure 
to  be  followed  is  given  in  Spicer's  report  C?J. 

6.  ACKNOWLSICEMENT 

Thanks  are  due  to  3.J.  Clement  for  his  help  with  the  graphical 
presentations . 


J.H.  James  (PSO) 


11  - 


REFERENCES 


1.  FULLER ,  C.R.,  FAHY  F.J.,  Characteristics  of  wave  propagation  anc. 
energy  distribution  in  cylindrical  elastic  shells  filled  with 
fluid,  J.  Sound  Vib.,  81(4),  1982,  pages  5OI-518. 

2.  JAMES,  J.H.,  Sound  radiation  from  fluid-filled  pipes,  Admiralty 
Marine  Technology  Establishment,  Teddington,  AMTE(N)TM81048, 

September  1981. 

3.  JAMES,  J.H.,  Computation  of  acoustic  power,  vibration  response 
and  acoustic  pressures  of  fluid-filled  pipes,  Admiralty  Marine 
Technplogy  Establisment,  Teddington,  AMTE(N)TM82036,  May  1982. 

4.  FAHY,  F.J.,  A  technique  for  measuring  sound  intensity  with  a 
sound  level  meter,  Noise  Control  Engineering,  9(3).  1977, 
pages  155-162. 

5.  FAHY,  F.J.,  Measurements  with  an  intensity  meter  of  the  acoustic 
power  jf  a  small  machine  in  a  room,  J.  Sound  Vib.,  57(3),  1978, 
pages  311-322. 

6.  REINHART,  T.E.,  CROCKER,  M.J.,  Source  identification  on  a  diesel 
engine  using  acoustic  intensity  measurements,  Noise  Control 
Engineering,  18(3),  1982,  pages  84-92. 

7.  SPICER,  W.J.,  Acoustic  intensity  vectors  from  an  infinite  plate 
with  line  attachments.  Admiralty  Marine  Technology  Establishment, 
Teddington,  AMTE(N)TM8l086,  October  1981. 

3.  LEISSA,  A.W.,  Vibration  of  shells,  NASA  SP-288,  1973. 

9.  SCHULTZ,  T.J.,  et  al,  Measurement  of  acoustic  intensity  in  reactive 
sound  field,  J.  Acoust.  Soc.  Am.,  57(6)  Part  1,  1975,  pages  1263-1265. 


-  13  - 


APPENDIX 


The  Spectral  Sis-placements 

The  excitation  is  either  a  radial  point  force,  of  magnitude  ?  , 
that  is  located  at  the  cylindrical  coordinates  'a.,2,2' ,  or  it  is  an 
interior  point  source  of  sound,  of  free-field  pressure  ?  exp(ik,R  /H3. 
that  is  located  at  (x  ,0,0). _  Xhe  spectral  iisplacements3are  obtained, 
via  Novozhilov  shell  theory  L3j,  as  the  solution  of_the  matrix  equation 
that  relates  displacements  to  external  excitations  [_2] 
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An  axial  point  force  excitation  is  obtained  by  __  setting  the  right-hand 
side  of  the  matrix  equation  to  w7oe^/2TTa,0,0_  . 
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FIG. 2  AXIAL  WAVENUMBER  VERSUS  FREQUENCY. 
SHELL  IN  VACUO  AAAA 
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FIG. 3  AXIAL  WAVENUMBER  VERSUS  FREQUENCY. 
SHELL  WITH  INTERIOR  WATER  EXTERIOR  VACUUM. 
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